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ABSTRACT 

We study chaos synchronization between Josephson junctions governed by a central junction with 
a time delay and demonstrate with numerical simulations the possibility of high quality synchro¬ 
nization. The results are important for obtaining high power systems of Josephson junctions, 
which are promising for practical applications. 
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I.INTRODUCTION 

Josephson junctions are paradigmatic nonlinear dynamical systems, see e.g. references in [1-2]. 
The study of chaos and its control in such systems [3-6] could be of huge practical importance from 
the application point of view. Chaos control in such systems is important for Josephson junction 
devices such as voltage standards, detectors, SQUID (Superconducting Quantum Interference De¬ 
vice), etc. where chaotic instabilities are not desired. Chaotic Josephson junctions [7-8] can be 
used for short-distance secure communications in wi-£ systems, high-bandwidth telecommunica¬ 
tions, ranging purposes, etc [9-10]. In a chaos-based secure communications a message is masked 
in the broadband chaotic output of the transmitter and synchronization between the transmitter 
and receiver systems is used to recover a transmitted message [3-6, 11-16]. 

It is well-known that Josephson junctions are a source of radiation with frequencies up to the 
terahertz (THz) region [17-18]. Terahertz waves have a wide range of possible uses, including 
security scanning, remote sensing chemical signatures of explosives, non-invasive applications in 
medicine [9]. One of the principal challenges in the held is to develop compact, low-cost, efficient 
THz sources. Additionally, for many remote sensing and imaging applications, high optical pow¬ 
ers are desirable, in part owing to the signihcant attenuation of THz radiation by water vapor in 
the atmosphere [19, 17-18]. Unfortunately the radiation from a single Josephson junction is very 
weak, of the order of pW. Synchronization of arrays of Josephson junctions is one of the ways to 
increase the radiation power from such sources [17, 19-20]. Upon achieving synchronization the 
radiation power will be proportional to the number of Josephson junctions squared. By synchro¬ 
nizing tens of thousands junctions radiation power as high as mW can be achieved. Recently it 
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was established that certain highly anisotropic cuprate superconductors naturally realize a stack 
of strongly electromagnetically coupled intrinsic Josephson junctions [17]. 

As shown in [21] the effect of the coupling delay time between junctions coupled in series on the 
synchronization quality is signihcant. As demonstrated by the numerical simulations in [21], with 
increasing coupling delays between the Josephson junctions correlation coefficients decay rapidly. 
The same tendency also holds for the increased number of junctions. However in [21] we were 
able to show that high quality synchronization is still a possibility among unidirectionally coupled 
in a series of 9 Josephson junctions on condition that the coupling delay between the junctions 
are reduced signihcantly (hundreds of times). It is clear that there is a physical limit to such a 
reduction. In light of this it is important to search for alternative ways of coupling between the 
junctions to synchronize very large number of the coupled systems. Results of numerical modelling 
presented in this paper show that coupling via a common medium-a central junction- could be 
one such viable alternatives. 

This paper is organized as follows. In section II we present a model of Josephson junctions 
coupled via a central junction with time delay. Section III deals with numerical simulations of the 
model under study and compares the result of simulations with the case of Josephson junctions 
coupled in series. Conclusions are presented in Section IV. Acknowledgements are presented in 
Section V. 


II. JOSEPHSON JUNCTIONS GOVERNED BY A CENTRAL JUNCTION 


Consider the following Josephson junctions 2,3,4 governed with a time delay by a central junction 
1 (Figure 1): 
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Where 0i, 02, 03, and 04 are the phase differences of the superconducting order parameter across 
the junctions 1,2,3 and 4, respectively;/! is called the damping parameter {j3R)^ = h{2eIcC)~^, 
where U, R and C are the junctions’ critical current, the junction resistance, and capacitance, 
respectively;/! is related to McCumber parameter /!c:/!^/!c = 1; is Planck’s constant divided by 
27r;e is the electronic charge;Zrfc is the driving the junctions direct current; It is worth noting that in 
general the total current passing through the junction will contain contributions from supercurrent 
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carried by the electron pairs, from a displacement current and an ordinary current. Representing 
the displacement current by a capacitor, and the ordinary current by a resistor one arrives at the 
equivalent circuit for Josephson junction shown in Figure 2. This junction circuit is a building 
block for the schematic circuits shown in Figure 1 and Figure 3. 

io cos(Qt + 0) is the driving ac (or rf) current with amplitudes fo,frequencies hi and phases 0;t is 
the coupling delay time between the junctions 1 — 2, 1 — 3 and 1 — 4; coupling between the junctions 
1 — 2,1 — 3 and 1 — 4 is due to the currents flowing through the coupling resistors Rs between the 
junctions 1 — 2,1 — 3 and 1 — 4 ;as = R/3Rj^ is the coupling strength between the junctions 1-2, 
1-3 and 1-4. We note that in Eqs.( 1-4)direct current and ac current amplitudes are normalized 
with respect to the critical currents for the relative Josephson junctions; ac current frequencies hi 
are normalized with respect to the Josephson junction plasma frequency = 2eIc{hC)~^, and 
dimensionless time is normalized to the inverse plasma frequency. 

We emphasize that Josephson junction 1 is the central one (see Figure 1) and is connected 
to junctions 2,3,4 unidirectionally through the coupling resistors Rg with some delay. We also 
underline there is no direct connection between the junctions 2,3,4 to be synchronized. This is 
a principal distinction from the connection topology considered in [21]. In [21] the Josephson 
junctions to be synchronized were connected unidirectionally in series (Figure 3). In Figures 1 
and 3 unidirectionality of the coupling between the junctions is provided by a Ferrite Isolator-FI. 
As the information propagation speed is hnite, any transmission (connection) line between the 
junctions will provide a delay time. 

It is noted that if the external drive current is purely dc then we have a second order au¬ 
tonomous system, Eq.(l) and therefore chaotic dynamic for Eq. (1) is ruled out. In order to make 
the dynamics of the Josephson junction (1) chaotic, one can add an ac current to the external 
driving, which makes Eq.(l) a second order non-autonomous system. Then treating Qt term as 
a new dynamical variable Eq.(l) (and consequently Eqs.(2-4)) can be rewritten as a system of 
third-order ordinary differential equations, which can behave chaotically. 
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It is worth noting that Eq.(l) is also used for modeling of a driven nonlinear pendulum, charge 
density waves with a torque and sinusoidal driving held [22-23]. 


III.NUMERICAL SIMULATIONS 

First we demonstrate that three Josephson junctions governed with time delay by the central 
junction can be synchronized. We note that the total number of junctions is four: a central 
junction connected to the other three remaining junctions (Figure 1). This is done in order to 
compare the results of this paper with the result of numerical simulations of synchronization 
between three connected in series junctions (Figure 3). We study chaos synchronization between 
Eqs.(5-16) using the correlation coefficient C [24] 

^ < i.x(t)- < X >){y{t)- <y>)> 

{x{t)- < X >)2 >< {y{t)- <y>f> 

where x and y are the outputs of the interacting systems; the brackets< . > represent the time 
average; this coefficient indicates the quality of synchronization: C = 1 means perfect synchro¬ 
nization. 

We simulate Eqs.(5-16) for the following set of parameters:/? = 0.25, = 0.3,0 = 0.7, = 

0.6,6* = 0,0;* = 0.45, r = 0.2. Parameter values used in the numerical simulations are taken 
from [19-20]. Figure 4 demonstrates the error dynamics between the junctions 2 — 4 (governed 
by junction 1) 04 — 02- While C 24 indicates the correlation coefficient between 02 and 04. These 
simulation results underline the high quality synchronization between three Josephson junctions 
2 — 4 with different initial conditions. Figure 4a shows chaotic dynamics of the central junction 
01 ( 0 . 

For comparison we also reproduce results of numerical simulations from [21] demonstrating 
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high quality synchronization between the end junctions coupled in series with a time delay. Fig¬ 
ure 5 shows the error dynamics between Josephson junctions 1 and 3 for the following 

set of parameters:/^ = 0.25, = 0.3, Zq = 0.7, hi = 0 . 6,6 = 0 ,as = 0.45, r = 0 . 2.(713 indicates the 

correlation coefficient between z/Zi and 

As mentioned in Section 1 the effect of the coupling delay time between the junctions cou¬ 
pled in series on the synchronization quality is significant: with increasing coupling delays be¬ 
tween the Josephson junctions correlation coefficients decay rapidly. The same tendency also 
holds for the increased number of junctions. However in [21] we were able to show that high 
quality synchronization is still a possibility among unidirectionally coupled in a series of 9 Joseph¬ 
son junctions on condition that the coupling delay between the junctions are reduced signif¬ 
icantly from r = 0.2 to r = 0.001. Figure 6 pictures z/Zg vs. -01 for the set of parameters 
/3 = 0.25, Zrfc = 0.3, Zq = 0.7, H = 0.6,9 = OjCr^ = 0.45, r = 0 . 001 .(Fig is the correlation coefficient 
between variables /Jg and z/Zi. 

Clearly there is a physical limit to such a reduction. If we keep r = 0.2 then the synchro¬ 
nization quality between the end junctions in an array of 9 junctions is relatively low. Figure 7 
presents the phase portrait between the dynamical variables /Jg and ijji for j3 = 0.25, idc = 0.3, zg = 
0.7,12 = 0.6,6 = 0,a:s = 0.45, r = 0.2.(Fig is the correlation coefficient between variables z/Zg and 
-01 • 

In the light of this, an alternative way of coupling between the junctions in order to synchro¬ 
nize a very large number of junctions investigated in this paper is of considerable importance. In 
Figure 8 we present the results of numerical simulations in an array of 30 Josephson junctions 
coupled via a common medium-a central junction (in total 31 junctions) and show high quality 
synchronization is possible between the junctions with random initial conditions. For Figure 8 
the coupling delay times between the central and other junctions to be synchronized are equal to 
r = 0.2. The other parameter values are as in Figure 4. C 2-31 is the correlation coefficient between 
variables Z //31 and z// 2 - As pictured in Figure 9 high quality synchronization takes place even in the 
case of 100 junctions coupled with the central junction with time delays r = 0 . 2 ; other parameter 
values are as in Figure 4. The correlation coefficient (F2-101 between dynamic variables z/zioi and 
ip 2 is equal to 0.99. 

IV.CONCLUSIONS 

In summary, we have numerically studied chaos synchronization between time-delay coupled 
Josephson junctions governed by a central junction. We have demonstrated the possibility of 
high quality synchronization between snch systems. We have shown that the topology of the 
connection between the jnnctions is a very important factor for the synchronization qnality. For 
many application pnrposes, e.g. in secnrity scanning, remote sensing of explosives, non-invasive 
medicine, high power Josephson jnnctions are reqnired. Synchronization of arrays of Josephson 
jnnctions is one way to increase the radiation power from snch sonrces. 
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FIGURE CAPTIONS 


FIG.1.Schematic sketch of three Josephson junctions 2,3,4 governed with time delay by the cen¬ 
tral junction 1 . Couplings 1 - 2 , 1-3 and 1-4 occur with the same strength as via a resistor Rg. 
Unidirectionality of couplings is provided by a Ferric Isolator (FI). 

FIG.2.Schematic set-up of the RC-shunted Josephson Junction subject to the external driving 
lext- The total current passing through the junction contains contributions from the supercurrent 
Is carried by electron pairs, the current through the capacitor C and the current through the 
resistor R. 

FIG.3.Schematic sketch of three Josephson junctions 1 — 2 — 3 coupled in series. Unidirectional 
couplings 1-2 and 2-3 occur with the same strength via a resistor Rg. Ferric Isolator (FI) 
provides unidirectionality of couplings. 

FIG.4.Numerical simulation of Eqs. (5-16): error dynamics of -04 — 'ip 2 '-C 2 A is the correlation coef- 
hcient between Josephson junctions 2 and 4. The parameters are:/I = 0.25, idc = 0.3, zq = 0.7, G = 
0.6, 0 = 0, = 0.45, r = 0.2.Dimensionless units. 

FIG.4a.Numerical simulation of Eqs.(5-7): Ghaotic dynamics of the central junction '0i(f). The 
parameters are:/5 = 0.25, Gc = 0.3, Zq = 0.7, D = 0.6, 6 * = 0, = 0.45.Dimensionless units. 

EIG.5.Numerical simulation of three Josephson junctions 1 — 2 — 3 coupled in series: error dy¬ 
namics of is the correlation coefficient between the end Josephson junctions 1 and 3. 

The parameters are:/? = 0.25, z^c = 0.3, Zq = 0.7, D = 0.6,6* = 6 , 0:5 = 0.45, r = 0.2.Dimensionless 
units. 

EIG. 6 .Numerical simulation of nine Josephson junctions coupled in series: -i/Zg vs.'i//i:Gi 9 is the 
correlation coefficient between the end Josephson junctions 1 and 9. The parameters are: /? = 
0.25, Zrfc = 0.3, Zq = 0.7, D = 0 . 6 , 6 * = 6 , 0:5 = 0.45, r = 0 . 001 .Dimensionless units. 

FIG.7.Numerical simulation of nine Josephson junctions coupled in series: Phase portrait between 
variables z/Zg andz/;i:Gi 9 is the correlation coefficient between the end Josephson junctions 1 and 9. 
The parameters are:/? = 0.25, z^c = 0.3, Zg = 0.7, D = 0.6,6 = 6 , 0:5 = 0.45, r = 0.2.Dimensionless 
nnits. 

FIG. 8 .Numerical simulation of 30 Josephson junctions governed by a central jnnction: error dy¬ 
namics of z/Zgi — 'il^ 2 '-C 2 - 3 i is the correlation coefficient between Josephson junctions 2 and 31. The 
parameters are: /? = 0.25, z^c = 0.3, Zg = 0.7, D = 0.6,6* = 0,o; = 0.45, r = 0.2. Dimensionless 
nnits. 

FIG.9.Numerical simulation of 100 Josephson jnnctions governed by a central junction: error dy¬ 
namics of z/Zioi — 'ip 2 '-C 2 -ioi is the correlation coefficient between Josephson jnnctions 2 and 101. 
The parameters are: /? = 0.25, idc = 0.3, zg = 0.7, D = 0.6, 6 = 0,as = 0.45, r = 0.2.Dimensionless 
units. 
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